Kolmogorov scaling and intermittency in Ray leigh- Taylor turbulence 
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The Rayleigh-Taylor (RT) turbulence is investigated by means of high resolution numerical sim- 
ulations. The main question addressed here is on whether RT phenomenology can be considered 
as a manifestation of universality of Navier-Stokes equations with respect to forcing mechanisms. 
At a theoretical level the situation is far from being firmly established and, indeed, contrasting 
predictions have been formulated. Our first aim here is to clarify the above controversy through a 
deep analysis of scaling behavior of relevant statistical observables. The effects of intermittency on 
the mean field scaling predictions is also discussed. 



The Rayleigh-Taylor (RT) turbulence is a well-known 
buoyancy induced fluid-mixing mechanism occurring in 
a variety of situations ranging from geophysics (see. e.g., 
Ref. [l[ in relation to cloud formation) to astrophysics (in 
relation to thermonuclear reactions in type-la supernovae 
0, Q and heating of solar coronal Q) to technological 
related problems, e.g., inertial confinement fusion (see 
Ref. i). 

Despite the ubiquitous nature of RT turbulence, a con- 
sistent phcnomcnological theory has been proposed only 
recently @. In three-dimensions, this theory predicts 
a Kolmogorov-Obukhov turbulent cascade in which tem- 
perature fluctuations are passively transported. This sce- 
nario, which is partially supported by numerical simula- 
tions [3j, |7j , has however been contrasted by an alternative 
picture which rules out Kolmogorov phenomenology Q . 

The goal of our work is twofold. From one hand we 
give stronger numerical support to the phenomenological 
theory a la Kolmogorov in RT turbulence. On the other 
hand, we push the analogy with usual Navier-Stokes (NS) 
turbulence much further: we find that small scale veloc- 
ity fluctuations in RT turbulence develop intermittent 
statistics analogous to NS turbulence. 

We consider the 3D, incompressible (V • v = 0), mis- 
cible Rayleigh-Taylor flow in the Boussinesq approxima- 
tion 
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where T is the temperature field, proportional to density 
via the thermal expansion coefficient f3, v the kinematic 
viscosity, k the molecular diffusivity and g = (0, 0, g) is 
the gravitational acceleration. 

At time t = the system is at rest with cooler (heavier) 
fluid placed above the hotter (lighter) one. This corre- 
sponds to v(x, 0) = and to a step function for the ini- 
tial temperature profile: T(x, 0) — — (6'o/2)sgn(z) where 
9q is the initial temperature jump which fixes the At- 
wood number A = (l/2)/3#o- The development of the 
instability leads to a mixing zone of width h which starts 



FIG. 1: (a) Snapshot of temperature field for Rayleigh-Taylor 
turbulence at t/r = 2.6. White (black) regions correspond to 
hot (cold) fluid, (b) Mean temperature profiles (T(z)) for 
times t/r = 1.4 (continuous), t/r = 2.0 (dashed) and t/r = 
2.6 (dotted), (c) Growth of the mixing layer thickness h(t) 
defined as the vertical range for which |(T(z))| < 0.98#o/2 
compensated with the dimensional prediction Agt 2 in order 
to get the dimensionless coefficient a. Filled symbols: a = 
h/{Agt 2 ), open symbols: a = h 2 /(AAgh) @|. 



from the plane z = and is dimensionally expected to 
grow in time according to h(t) = a Agt 2 [1, [§]. Inside 
this mixing zone, turbulence develops in space and time. 
The phenomenological theory Q predicts for velocity and 
temperature fluctuations the scaling laws 



5 r v(t) ~ (Agflh^r 1 / 3 
5 r T(t) ~ OoiAgyVh-Wr 1 ' 3 



(3) 
(4) 



The first relation represents Kolmogorov scaling with a 
time dependent energy flux e ~ (Ag) 2 t. ^From the above 
scaling laws one obtains that the buoyancy term (3gT 
becomes subleading at small scales in ([T]), consistently 
with the assumption of passive transport of temperature 
fluctuations. 
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We integrate equations ([T][2]) by a standard 2/3- 
dealiased pseudospectral method on a periodic domain 
with uniform grid spacing, square basis L x = L y and as- 
pect ratio L x /L z — r with a resolution up to 512 x 512 x 
2048 (r — 1 /4) . Time evolution is obtained by a second- 
order Runge-Kutta scheme with explicit linear part. In 
all runs, Ag = 0.25, Pr = v/n = 1, 0q = 1. Viscosity is 
sufficiently large to resolve small scales (k max r\ ~ 1.2 at 
final time). In the results, scales and times are made di- 
mensionless with the box scale L z and the characteristic 
time r = {L z /Ag) 1 / 2 0. 

Rayleigh- Taylor instability is seeded by perturbing the 
initial condition with respect to the step profile. Two dif- 
ferent perturbations were implemented in order to check 
the independence of the turbulent state from initial con- 
ditions. In the first case the interface T = is perturbed 
by a superposition of small amplitude waves in a narrow 
range of wavenumber around the most unstable linear 
mode [n|. For the second set of simulations, we per- 
turbed the initial condition by "diffusing" the interface 
around z — 0. Specifically, we added a 10% of white 
noise to the value of T(x, 0) in a small layer of width ho 
around z = 0. 

Figure [T] shows a snapshot of the temperature field for 
a simulation with r = 1/2 at advanced time. Large scale 
structures (plumes) identify the direction of gravity and 
break the isotropy. Nonetheless, we find that at small 
scales isotropy is almost completely recovered: the ra- 
tio of vertical to horizontal rms velocity is v z /v x ~ 1.8 
while for the gradients we have d z v z /d x v x ~ 1.0. The 
horizontally averaged temperature (T(z)) follows closely 
a linear profile within the mixing layer where, therefore, 
the system recovers statistical homogeneity. 

The analysis of the mixing layer width growth is also 
presented in Fig. [T] As shown by previous studies d, [l2| , 
the naive compensation of h(t) with Agt 2 does not give 
a precise estimation of the coefficient a because of the 
presence of subleading terms which decay slowly in time. 
We have therefore implemented the similarity method 
introduced in Q which gives an almost constant value of 
a ~ 0.038 for t/r > 1.5, consistent with previous studies 

Figure [2] shows the kinetic energy E(k) and tempera- 
ture Ex{k) spectra within the similarity regime. ^From 
([3]) and (|U), we expect the following spatial-temporal 
scaling of spectra: E(k,t) ~ t 2 / 3 fc~ 5 / 3 and E T {k,t) ~ 
£ -4 / 3 fc -5 / 3 . Kolmogorov scaling fc~ 5 / 3 is evident for both 
velocity and temperature fluctuations. Moreover, self- 
similar temporal evolution of spectra is well reproduced, 
as shown in the lower inset. Also in Fig. [5] the two con- 
tributions to kinetic energy flux in spectral space are 
shown. Buoyancy contribution, dominant at large scale, 
becomes subleading at smaller scales, in agreement with 
the Kolmogorov-Obukhov picture. The above results, 
together with previous simulations 0, @] and theoretical 
arguments [y], give a coherent picture of RT turbulence 
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FIG. 2: Two-dimensional kinetic energy spectrum (o) and 
temperature spectrum (A) at time t/r = 2.6 corresponding 
to R\ — 245. Spectra are computed by Fourier transforming 
velocity and temperature fields on two-dimensional horizon- 
tal planes and averaging over z in the mixing layer. Dashed 
lines represent Kolmogorov scaling fc -5 ^ 3 . Lower inset: evo- 
lution in time of the amplitude of kinetic energy (x) and 
temperature (+) spectra at fixed wavenumber fco = 12. Lines 
represent the dimensional predictions t 2 ^ 3 (continuous) and 
t~ A/z (dashed) given by <(3j and Q. Upper inset: inertial 
(continuous) and buoyancy (dashed) contributions to kinetic 
energy flux Il(fc) in Fourier space. 



as a Kolmogorov cascade of kinetic energy forced by large 
scale temperature instability. 

In the following we push this analogy one step ahead 
by showing that small scale fluctuations in RT turbulence 
display intermittency corrections typical of usual Navier- 
Stokes (NS) turbulence. Intermittency in turbulence is 
a consequence of non-uniform transfer of energy in the 
cascade which breaks down scale invariance. As a conse- 
quence, scaling exponents deviates from mean field the- 
ory and cannot be determined by dimensional arguments 
[13j. Several studies have been devoted to the intermit- 
tent statistics in NS turbulence, where the main issue 
concerns the possible universality of anomalous scaling 
exponent with respect to the forcing mechanisms and the 
large scale geometry of the flow. While universality has 
been demonstrated for the simpler problem of passive 
scalar transport, it is still an open issue for nonlinear 
NS turbulence. Therefore the key question is whether 
small scale statistics in RT turbulence is equivalent to 
the statistics observed in homogeneous, isotropic turbu- 
lence. 

The simplest, and historically first, evidence of inter- 
mittency is in the dependence of energy dissipation on 
Reynolds number [3, EH, [lj| . Classical statistical indi- 
cators are the flatness K of velocity derivatives HI, 3] 
(corresponding to K ~ (e 2 )/ (e) 2 in terms of energy dissi- 
pation), and the variance of the logarithm of kinetic en- 
ergy dissipation which is expected to grow with Reynolds 
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number as 



cine = a + (3/u/2) ]nRx 



(5) 



The exponent \x is the key ingredient for the log-normal 
model of intermittency and its value is determined exper- 
imentally 15|, ll7[ and numerically [l8| to be [i ~ 0.25. 
More in general, moments of local energy dissipation are 
expected to have a power-law dependence on R\ 



(0) 



where the set of exponents t p can be p redicted within the 

lE ES in terms of 



multifractal model of turbulence 
the set of fractal dimensions D(h). 
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FIG. 3: Scaling of the variance of hie on Reynolds number 
defined as R\ = {v z ) rma {d z v z ) rms /v, obtained from two real- 
izations with white noise initial perturbation. The line is the 
best fit corresponding to /j, = 0.24 in ([5]). Inset: scaling ex- 
ponents of the moments of local dissipation t p obtained from 
best fits according to {BJ. The line represents the log-normal 
approximation t p — (3/4)/i(p 2 — p). 

Because in RT turbulence the Reynolds number in- 
creases in time, it provides a natural framework for a 
check of ([5]) and © ■ Figure [3J shows the dependence of 
the variance of In e on R\ together with the first moments 
of energy dissipation. Despite the limited range of R\ , a 
clear scaling of In e is observable, even if with some fluctu- 
ations. The best fit with © gives an exponent [i ~ 0.24, 
very close to what observed in homogeneous, isotropic 
turbulence [18 1. 

Scaling exponents t p for the moments of dissipation 
^ are also shown in Fig. [3J We were able to com- 
pute moments up to p — 2 with statistical significance. 
Log-normal approximation, which is in general valid for 
p — > 0, is found to be unsatisfactory for larger values 
of p. For p = 2, which corresponds to the flatness K 
of velocity derivatives, we find T2 ~ 0.27. This result 
is consistent with experiments at comparable Reynolds 
numbers [HI which shows that K ~ R\ 2 for R\ < 200 
while an asymptotic exponent ti ~ 0.41 is reached for 
Rx > 10 3 only. 




FIG. 4: Scaling exponents of isotropic longitudinal velocity 
structure functions S p (r) — {(5 r v ■ f) p ) (f = r/r) for the late 
stage of RT turbulence (open circle). Exponents are com- 
puted by compensation of S p (r) with Ss(r), according to the 
ESS procedure [2l|] averaging inside the mixing layer and on 
all directions. Filled circles: scaling exponents from simula- 
tions of homogeneous, isotropic turbulence at R\ = 381 [221 ], 
Line represents dimensional prediction £ p — p/3. Inset: third- 
order isotropic longitudinal structure function Ss(r). The line 
represents Kolmogorov's four-fifth law Ss(r) — — 4/5er. 



In NS turbulence intermittency is also observed in the 
inertial range of scales as deviations of velocity structure 
functions S p (r) — ((5 r v) p ) from the dimensional predic- 
tion j3j which corresponds to S p (r) ~ r p / 3 fl3j | ■ Anoma- 
lous scaling is observed, which corresponds to scaling 
laws S p (r) ~ with a set of exponents ( p ^ p/3. 
We remind that constancy of energy flux in the iner- 
tial range implies £3 — 1 independently on intermit- 
tency, as required by the Kolmogorov's "four-fifths" law 
S3 (r) = — 4/5er [l3| ]. which is indeed observed in our 
simulations (see inset of Fig. 0]) . Figure 2] shows the first 
longitudinal scaling exponents computed from our simu- 
lations exploiting the extended self-similarity procedure 
which allows for a precise determination of the exponents 
at moderate Reynolds numbers [2l| . A deviation from 
dimensional prediction £ p = p/3 is clearly observable for 
higher moments. Fig. [4] also shows the scaling exponents 
obtained from a homogeneous, isotropic simulation of NS 
equations at a comparable R\ The two sets agree 

within the error bars, this gives further quantitative evi- 
dence in favor of the equivalence between RT turbulence 
and NS turbulence in three dimensions. 

We end by discussing the behavior of turbulent heat 
flux and rms velocity fluctuations as a function of the 
mean temperature gradient. In terms of dimcnsionlcss 
variables, these quantities are represented respectively 
by the Nusselt number Nu = 1 + (v z T)L/(k9 ), the 
Reynolds numbers Re — v rms L/v and Rayleigh number 
Ra = AgL 3 / '(vk). The relations between these quanti- 
ties has been object of many experimental and numerical 
studies in past years, mainly in the context of Rayleigh- 
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Benard turbulent convection 23, 24, \2 



27, m 



Experiments have reported both simple scaling laws 
Nu ~ Rci 13 with exponent Q scattered around 3 = 0.3 
[29l | and, more complicated behavior [2^, [3(| par- 
tially in agreement with a phenomenological theory 24J. 
However, in the limit of very large Ra, Kraichnan [3l| 
predicted an asymptotic scaling Nu ~ Ra 1 / 2 now called 
the ultimate state of thermal convection. This regime is 
expected to hold when thermal and kinetic boundary lay- 
ers become irrelevant, and indeed has been observed in 
numerical simulation of thermal convection at moderate 
Ra when boundaries are artificially removed [27]. It is 



which turns out to be an excellent setup for experimental 
studies in this direction. 



therefore natural to expect that the ultimate state scal- 
ing arises in RT convection where boundaries are absent. 

The ultimate state relations can formally be obtained 
from kinetic energy and temperature balance equations 
24J. In the context of RT turbulence, they are a sim- 
ple consequence of the dimensional scaling of the mixing 
length L = h ~ Agt 2 and of the rms velocity v rms ~ Agt. 
Inserting in the definition of the dimensionless numbers 
one obtains 

Nu ~ Pr 1/2 Ra 1/2 , Re ~ P r - 1/2 Ra 1/2 (7) 

where Pr = u/k. 
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FIG. 5: The scaling of Nusselt number (open circles) and 
Reynolds number (solid circles) as functions of Rayleigh num- 
ber. Lines represent the ultimate state predictions @. 

We remark that the above relations are independent on 
the statistics of the inertial range and on the presence of 
intermittency as they involve large scale quantities only. 
Our numerical results, shown in Fig. [5l confirms the ul- 
timate state scaling ([7]). The same behavior has been 
predicted and observed for two-dimensional RT simula- 
tions, where temperature fluctuations are not passive and 
Bolgiano scaling is observed in the inertial range [jjj ]. 
The elusive Kraichnan scaling in thermal convection finds 
its natural manifestation in Rayleigh- Taylor turbulence, 
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